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Abstract: The performance of implantable electrochemical glucose sensors is highly 
dependent on the flux-limiting (glucose, H2O2, O2) properties of their outer membranes. A 
careful understanding of the diffusion profiles of the participating species throughout the 
sensor architecture (enzyme and membrane layer) plays a crucial role in designing a robust 
sensor for both in vitro and in vivo operation. This paper reports the results from the 
mathematical modeling of Clark's first generation amperometric glucose sensor coated 
with layer-by-layer assembled outer membranes in order to obtain and compare the 
diffusion profiles of various participating species and their effect on sensor performance. 
Devices coated with highly glucose permeable (HAs/Fe ) membranes were compared with 
devices coated with PSS/PDDA membranes, which have an order of magnitude lower 
permeability. The simulation showed that the low glucose permeable membrane 
(PSS/PDDA) sensors exhibited a 27% higher amperometric response than the high glucose 

">_|_ 

permeable (HAs/Fe ) sensors. Upon closer inspection of H2O2 diffusion profiles, this 
non-typical higher response from PSS/PDDA is not due to either a larger glucose flux or 
comparatively larger O2 concentrations within the sensor geometry, but rather is attributed 
to a 48% higher H2O2 concentration in the glucose oxidase enzyme layer of PSS/PDDA 
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coated sensors as compared to HAs/Fe coated ones. These simulated results corroborate 
our experimental findings reported previously. The high concentration of H2O2 in the 
PSS/PDDA coated sensors is due to the low permeability of H 2 0 2 through the PSS/PDDA 
membrane, which also led to an undesired increase in sensor response time. Additionally, it 
was found that this phenomenon occurs for all enzyme thicknesses investigated (15, 20 and 
25 nm), signifying the need for a holistic approach in designing outer membranes for 
amperometric biosensors. 

Keywords: amperometric glucose sensors; mathematical modeling; biosensors; 
Michaelis-Menten constant; enzyme kinetics; layer-by-layer assembly; polymer chemistry; 
finite difference schemes; differential equations; analyte diffusion 



1. Introduction 

The design and implementation of miniaturized implantable potentiometric or amperometric 
sensors for glucose monitoring has recently gained considerable attention in light of their ability to be 
integrated with integrated electronic circuits in the development of continuous glucose monitoring 
systems [1-6]. Among these, Clark's type amperometric sensors drew a lot of attention which is owed 
to their design simplicity and ease of fabrication. First generation Clark-based amperometric sensors 
employ a glucose oxidase (GO x ) enzyme immobilized on top of the working electrode [7]. The flavin 
adenine dinucleotide (FAD) redox cofactor of GO x catalyzes the oxidation of glucose to 
glucarolactone, as shown in Equations (1) and (2): 

Glucose + GO x (FAD) -> Glucolactone + GO x {FADH 2 ) (1) 

GO x {FADH 2 ) + 0 2 ^ GO x {FAD) + H 2 0 2 . (2) 

The generated H2O2 is amperometrically assessed on the surface of working electrode by the 
application of the appropriate redox potential, hence relating the current to glucose concentration, 

+v 

H 2 0 2 ^0 2 +2e~ + 2H+ (3) 

As evident from Equations (1) and (2), optimum sensor performance can only be obtained when the 
glucose to oxygen ratio is less than 1 . If not, Equation (2) becomes oxygen limited, which is the case 
for most in vivo applications [8]. For example, in the subcutaneous tissue the O2 concentration is only 
0.18 mM as compared to 5.6 mM of physiological glucose concentration (glucose/02 ratio = 30) [9], 
which leads to signal saturation at higher glucose concentrations. In order to mitigate signal saturation 
and render the sensor response linear, perm-selective membranes based on pore size has been 
employed, which decrease the glucose to O2 ratio in the vicinity of GO x (Figure 1). However, the use 
of outer membranes inadvertently decreases sensitivity and increases response time. 

Second and third generation Clark-type biosensors employ redox mediators and direct 'wiring' of 
enzymes to electrodes in an attempt to minimize or completely eliminate the need for O2 [9]. However, 
the unwanted toxicity associated with the nanomaterials used in second and third generation Clark type 
glucose sensors stand against their application in implantable sensor configurations [9]. With this in 
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mind, our group has been evaluating mediator-free, 1st generation Clark-type sensors that solely rely 
on careful mass transfer balance of the various species involved in the operation of these devices [10-13]. 
As part of this, using a series of layer-by-layer assembled glucose flux limiting membranes (namely 
humic acids/ferric cations (HAs/Fe 3+ ) and poly(styrene sulfonate)/poly(diallyldimethylammonium chloride) 
(PSS/PDDA)), we have shown that the outer diffusion of H 2 C>2 through these membranes plays a vital 
role in determining the sensitivity of the sensor. It was proven that sensors fabricated with highly 
glucose permeable HAs/Fe 3+ membranes exhibited lower sensitivities and better linearities compared 
to sensors with less permeable PSS/PDDA membranes [10]. Through a series of analyte diffusion 
studies and oxygen dependence of these sensors [10], it was proven that this unique phenomenon did 
not occur as a result of the poor diffusion of oxygen through the outer membranes, but rather as a 
result of the outer diffusion of H 2 C>2 through the various outer layers. In light of these findings, we 
have shown that optimal sensor performance could be easily achieved via a careful stratification of 
five functional layers, each of which is intended to carefully balance the mass transfer of various 
participating species involved [11]. This has also prompted us to theoretically understand the diffusion 
profiles of various participating species as they diffuse through the polymer and enzyme layers, which 
is a powerful tool useful in facilitating the development of future high performance first generation 
Clark-based amperometric implantable sensors. 



Figure 1. Cross-sectional schematic representation of the multi-layer amperometric 
glucose sensor used in the numerical simulations. The outer LBL membrane is represented 
as either one of the following systems: (i) HAs/Fe (ii) PSS/PDDA. 
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Herein, we present our initial results from our theoretical analysis of glucose sensors employing the 
two aforementioned LBL membranes. A multi-layer analyte diffusion model based on Fick's second 
law was employed to analyze the diffusion profiles of the participating species (namely, glucose and 
H2O2). The computational model presented here utilizes well-developed enzyme-catalyzed reaction 
diffusion equations, which were applied to glucose and H2O2, to model the reaction of glucose with 
glucose oxidase. This along with Fick's second law of diffusion, which represents the diffusive 
transport of glucose and the enzyme-generated H2O2 throughout the LBL assembled outer membranes, 
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gives a full dynamic representation of the sensor-system under study. The explicit finite difference 
scheme was coded in MATLAB to solve the diffusion-based partial differential equations in order to 
obtain concentration profiles of glucose and H2O2 as a function of LBL thickness and enzyme 
thickness. In doing this, we were able to precisely analyze the role of H 2 0 2 in sensor operation and 
its effect on the amperometric response. In accordance with our previously reported experimental 
studies [10], sensors coated with the less glucose permeable PSS/PDDA membrane showed a 
pronounced build-up of H2O2 compared to that of sensors coated with the more permeable HAs/Fe 
membrane. It was also found that this phenomenon appears for increasing enzyme thicknesses, further 
confirming the previous experimental studies. Furthermore, amperometric response times were found 
to increase with increasing enzyme thickness, due to the time necessary for the participating species to 
establish equilibrium within the enzyme layer. 

2. The Theoretical Model 

Figure 1 shows the cross-section of the sensor under study. The platinum electrode is covered with 
an electropolymerized GO x /poly(o-phenylenediamine) (PPD) layer followed by a layer-by-layer 
assembled outer membrane based on HAs/Fe (Device A) and PSS/PDDA (Device B). As can be seen 
in Figure 1 , the incorporation of GO x into the electropolymerized poly(o-phenylenediamine) (PPD) has 
been previously shown to block the oxidation of other endogenous species like uric acid, ascorbic acid 
and acetaminophen, which are highly likely to oxidize at the operating potential of the sensor [14-16]. 
The dynamic operation of this multi-layer sensor can be accurately modeled by one-substrate 
Michaelis-Menten enzyme kinetics [17-21]. First, consider the enzyme-catalyzed reaction where the 
substrate, Se, is irreversibly converted to a product, P E : 

ki k2 

S E + E <± ES E -> E + P E (4) 
k-i 

Here, the substrate, S E , is the glucose present in the enzyme layer, the product, P E , is the hydrogen 
peroxide (H2O2) in the enzyme layer generated per Equation (2), and ki, k-i and k 2 are the rate 
constants. The rate at which the product is formed, V, can be expressed as: 

V = ^ = k 2 [ES E ] (5) 

Applying the quasi steady-state assumption on the rate of variation of the [ESe] complex [17] leads 
to the following equation: 

v = g[Pg] = v m [Se] (6) 

dt K M +[S E ] ^ ' 

where Km is the Michaelis-Menten constant and Vm is the maximum enzymatic rate when the enzyme 
is fully saturated. We can integrate this rate equation to the well-known diffusion equation to obtain 
expressions for the effect of enzyme-catalyzed reactions on both glucose and H2O2 as they diffuse 
through the GO x /PPD layer (layer 1). Here we can safely assume that the GO x is uniformly distributed 
throughout the layer [22], which is a self-limiting electropolymerization process that yields a thickness 
comparable to the size of the immobilized enzyme [23]. Additionally, we assume that oxygen is 



Sensors 2012, 12 



13406 



present in sufficient amounts to drive the enzematic reaction to completion. The results in two 
reaction-diffusion equations (Equations (7) and (8)): 



dS E 
dt 

dP E 
dt 



D 



d 2 S E 



* S E 



S E d> 



,0 < x < d E ,0 <t <T 



= D 



P E dx 2 



E + Vjna^E iQ <X <d E ,0 <t<T 



(V) 
(8) 



K M +s E 

where S E (x, t) and P E (x, t) are the concentrations of glucose and H2O2, respectively, D S e the diffusion 
coefficient of glucose in the GO x /PPD layer, Dp E the diffusion coefficient of the H2O2 in the GO x /PPD 
layer, and d E is the total thickness of the GO x /PPD layer. Here, x and t represent the points in space 
and time, respectively, where x = 0 is taken at the electrode surface, and x = d E at the edge of the 
GOx/PPD layer, and T is the total time elapsed. A graphical representation of the boundary conditions 
of this system conditions can be seen in Figure 2. 

Figure 2. Graphical representation of the Product (H2O2), P, and Substrate (Glucose), S, 
boundary conditions employed in the simulated multi-layered glucose sensor. Here, x = 0 
is taken at the surface of the Pt electrode, x = is taken as the enzyme-membrane 
interface and x = ds +dM is taken as the surface of the LBL outer membrane. 
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The operation of the system initiates immediately after the substrate in the bulk solution reaches the 
surface of the outer membrane (layer 2), which we define as the following initial conditions: 

S B (x, 0) = 0, P E (x, 0) = 0, x G [0, d E ] (9) 

P M (x, 0) = Ox G [d E , d E + d M ] (10a) 

S M {dE, 0) = 0 S M {dE + dM, 0) = S 0 (10b) 

Here, we introduce S M (x,t) and P M (x,t) which are the concentrations of glucose and H 2 0 2 , 
respectively, in the outer LBL membrane, where x = d M represents the surface of the outer membrane. 
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The assumptions made are the following: (1) the substrate does not react at the electrode surface, 
(2) the H2O2 at the electrode surface gets permanently reduced to zero, (3) if the substrate is well 
stirred, then the diffusion layer (0 < x < d E + d M ) remains at a constant thickness, d E + d M , during 
sensor operation and (4) the concentration of H 2 0 2 at the sensor/bulk solution interface is zero. This 
can be quantitatively described as the following boundary conditions: 



dS E 
dx 



= 0 (11) 

S M (d E + d M ,t) = S 0 (12) 

P £ (0,t) =P M {d E + d M ,t) = 0 (13) 

The aforementioned initial conditions and boundary conditions have been previously employed to 
model various amperometric and potentiometric biosensors, with various physical parameters and 
testing conditions [18-21,24-31]. At the boundary between the enzyme layer and the flux-limiting 
membrane, matching conditions are defined (for t > 0) [32]: 



S E dx 



n Mb 

e ax 



— D SSm 

x=dE ~ S M dx 



,S E {d E ,t) = S M {d E ,t) (14) 

x=d E 



x=d E Pm 



,P E (d E ,i) = P M {d E ,t) (15) 

x=d E 



Following H 2 0 2 oxidation at the Pt surface, the amperometric current response of the device is 
expressed as: 

,dP E 



I = n e D PE F — 



(16) 

x=0 



where, n e is the number of electrons produced (2 for H 2 0 2 ), D PE the diffusion coefficient of the H 2 0 2 
in the GO x /PPD layer and F is Faraday's constant (96,500 C/mol). 

2.2 The Solution Scheme 

The cross-sectional schematic of Figure 2 was modeled by constructing a mesh employing a finite 
set of x and t, and utilizing the boundary conditions of Equations (11) to (15) as a means to obtain a 
solution to the tridiagonal system. Based on this, the explicit finite difference scheme was used to 
obtain substrate and product concentrations of the Michaelis-Menten reaction-based differential 
Equations (7) and (8), as: 

Glucose (S) 

^i +1 = 4+^(4 -24 + 4 )-At^4 

h i h i Ax 2 v h i+i h i t i-i J Km+S j 

E i 

H2O2 (P) 

Pi +1 =Pl+^{Pl -H+Pl )+At^4 (18) 

E i 

where 7 = 1 ... Numt, where Numt represents the total number of time steps and i = I ... k— I, where k 
represents the number of steps in space from the electrode surface to the enzyme/LBL interface. Having 
said this, at the interface (i = k) Equations (14) and (15) can be solved explicitly [32] as follows: 
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Glucose (S) 

= S ! + £- 2 WsmSL ~ (D SM + D SE )S( + D SE SU] (19) 

H 2 0 2 (P) 

P/ +1 = Pi + £ [Wtii - OW + D„ )?/ + Dp^P/lj (20) 

where i = k- I ... k+ I. Furthermore, similar to the difference scheme of Equations (17) and (18) and 
neglecting the Michaelis-Menten reaction portion, the diffusion of glucose and H 2 C>2 from the bulk 
through the LBL membrane follows the diffusion equation and was solved as,: 

Glucose (S) 

= sL, + S g £ GL-2S£, + iiJ (17) 

H20 2 (P) 

where i = k + 1 ... Numx, where Numx represents the total number of time steps performed in the 
simulation. 

Numerical parameters for the enzymatic layer, such as diffusion coefficients, Michaelis-Menten 
constant, maximum rate of velocity and thickness have been obtained from earlier reports [22,33]. A 
keen understanding of the kinetics of H2O2 diffusion though this multi-layer model has been simulated 
and analyzed by studying the effect of enzyme thickness (15, 20 and 25 nm GO x /PPD) as a function of 
these two LBL outer membranes. A summation of variables employed in the model can be seen in 
Table 1. 



Table 1. Compilation of variables used in the multi-layer model. 



Declaration of Variables 


S E 


Concentration of glucose in the 
GO x /PPD layer 


Dsm 


Diffusion coefficient of glucose in the various LBL 
outer membranes 


Pe 


Concentration of H 2 0 2 in the GO x /PPD 
layer 


Dpm 


Diffusion coefficient of H 2 0 2 in the various outer 
membranes 


Sm 


Concentration of glucose in the various 
LBL outer membranes 


Vmax 


Maximum enzymatic rate of reaction 


Pm 


Concentration of H2O2 in the various 
outer membranes 


Km 


Michaelis-Menten constant 


d E 


Thickness of the GO x /PPD layer 


I 


Response current 


d M 


Thickness of the various LBL outer 
membranes 


n e 


Number of electrons (2 for H 2 0 2 ) 


D SE 


Diffusion coefficient of glucose in the 
GO x /PPD layer 


F 


Faraday's Constant (96,500 C/mol) 


Dpe 


Diffusion coefficient of H 2 0 2 in the 
GO x /PPD layer 
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The physical parameters for the enzyme layer, as well the parameters of the two outer membranes 
are summarized in Table 2. 



Table 2. Physical parameters for the enzyme layer and LBL assembled outer membrane 
used in the simulation, including diffusion coefficients for glucose and H 2 0 2 , membrane 
thicknesses, Michaelis-Menten constant and maximum rate of enzymatic activity. 





Layer 


Parameter 


GO x /PPD 


HAs/Fe 3+ 


PSS/PDDA 


D G i ucose (crn/sec) 
D H20 2 (cm 2 /sec) 
Thickness (nm) 
K M (mol/cm 3 ) 
Vmox {mol/cm 2 sec) 


4 x 10" 10 [22,33] 

5 x 1(T 9 [22,33] 
15, 20 and 25 [33] 
33 x 10" 6 [22,33] 
60 x 10" 6 [22,33] 


4.4 x 10" 9 [10] 
4.8 x 1(T 9 [10] 
100 [10,12] 

n/a 
n/a 


0.058 x 10" 9 [10] 
0.1 x 10 9 [10] 
5 [10,12] 

n/a 
n/a 



3. Results 



It is well believed that the amperometric response of first generation Clark's glucose sensors is 
directly proportional to the glucose permeability through the outer membranes assuming that O2 is 
not the limiting factor for the GO x reaction (Equation (2)). Having said this, in order to keep the 
glucose-to-oxygen ratio less than 1, a variety of outer membranes with less and less glucose 
permeability have been reported [8,9], albeit without a mechanistic understanding of the diffusion 
profiles of other participating species such as H2O2. Layer-by-layer assembled outer membranes afford 
such studies in that the thickness as well as analyte permeability can be precisely controlled via the 
number of bilayers within the membrane. Utilizing a series of such LBL assembled membranes 
(namely a more glucose permeable HAs/Fe 3+ and less glucose permeable PSS/PDDA membrane), we 
have previously deduced that the amperometric response of a sensor is also dependent on the outward 
diffusion of H2O2. Intuitively, the more glucose-permeable HAs/Fe 3+ outer membrane is expected to 
yield a higher response; yet the outwards diffusion of H2O2 through these membranes is also 
significant thereby contributing to a lowering of the amperometric response. In addition, through 
a series of oxygen dependence studies it was concluded that any potential role of lower oxygen 
levels within HAs/Fe membranes in lowering the amperometric response of HAs/Fe device has 
been excluded. 

In this contribution, we report a simulation on the concentration profiles of both glucose and H2O2 
throughout the multi-layer system in an attempt to quantify their peak concentrations within the 
interior of the sensor and to re-assert our earlier finding that H2O2 is a key contributor in advanced 
implantable sensor design. Here, it is worth noting that the theory and mathematical analysis of 
enzyme-catalyzed reactions used in this model are widely-accepted, and thus have been previously 
applied to similar glucose oxidase systems along with a wide-range of other enzymatic processes. 
More importantly, based on our earlier experimental studies on oxygen dependence of implantable 
glucose sensors, we have in this study assumed that O2 is present in sufficient amounts (for 
Equation (2) to go to completion) and that there is no variation in oxygen levels during sensor 
operation. Moreover, any effect of oxygen within the electropolymerized PPD/GOx layer can be 
considered negligible based on the fact that the pore size of PPD is enough to let H2O2 pass through, 
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which is a bigger molecule than Oj. In fact, it has been reported that sensors having PPD/GO x 
structures are minimally affected by oxygen [34], reinforcing our 02-independnce assumption. Having 
said this, it should be clear that the assumption of abundant O2 levels within the sensor is an 
approximation, and has been used in order to solely analyze the role of H2O2 in these devices. 

Based on that, we have utilized a one-substrate Michaelis-Menten model [18] in conjunction with 
diffusion kinetics through the LBL assembled outer membranes to get a holistic representation of 
sensor performance for the comparison of both devices. 

Figure 3 shows the simulated amperometric response of sensors coated with the two LBL systems. 
For any glucose concentration, and in particular at higher glucose concentrations, the amperometric 
response of Device B (with PSS/PDDA membranes) is 27% higher than that of Device A (HAs/Fe 3+ 
membranes). This trend in amperometric response is opposite to the trends observed in glucose 
permeability through the LBL membranes (which is highest through the HAs/Fe , Device A, and least 
through the PSS/PDDA, Device B, membranes), wherein the amperometric response from Device A is 
27-30% lower than Device B. Here, it is worth noting that the trend in simulated amperometric 
responses is similar to that of experimental data reported earlier [10]. Taking into account the low 
oxygen permeability of PSS/PDDA membranes compared to HAs/Fe 3+ membranes, it was concluded 
that the observed higher amperometric response from Device B (PSS/PDDA) is due to a higher 
build-up of H2O2 within its sensor geometry. The higher build-up of H2O2 is in turn attributed to its 
low permeability through the PSS/PDDA membrane compared to the HAs/Fe membrane. In order to 
further elucidate this, a theoretical understanding of the concentration profiles of glucose and H2O2 as 
a function of various LBL membranes has been performed. 

Figure 3. Simulated amperometric response of Device A (100 nm HAs/Fe 3+ ) and Device B 
(5 nm PSS/PDDA) employing 20 nm GO x /PPD coating the electrode surface. 
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Figure 4(A,B) illustrates the simulated concentration profile of glucose and H2O2 within a sensor 
geometry coated with either of the two LBL membranes namely (HAs/Fe 3+ and PSS/PDDA) at normal 
physiological glucose concentrations (6 mM glucose). The GO x enzyme for both devices is 
immobilized within a 20 mM poly(phenylenediamine) (PPD) layer and is assumed to be uniformly 
distributed within the PPD layer [22]. Employing the outer membrane thickness values of Table 1, the 
bulk solution of Device A is present at 120 nm from the electrode surface, while the bulk solution is 
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present at 25 nm from the electrode surface for Device B. As illustrated in Figure 4(A,B), the glucose 
flux within the GO x layer for Device B (PSS/PDDA) lower than that in Device A (HAs/Fe 3+ ), which is 
in agreement with the trends in glucose permeability through these membranes [12]. Based on the 
observed lower glucose flux within Device B, one would expect it to have a lower concentration of 
H 2 0 2 generated compared to Device A. However, as seen in Figure 4, the concentration of H 2 0 2 
within Device B is 48% higher than that of Device A. Considering the facts that the rate of oxidation 
of H 2 0 2 at the Pt electrode is the same for all the electrodes and that the oxygen permeability through 
the three membranes follows the same trend as glucose [10], the higher concentration of H 2 0 2 within 
Device B can easily be attributed to the poor outwards diffusion of H 2 0 2 through the PSS/PDDA 
membrane compared to the HAs/Fe 3+ membrane. This larger build-up of H 2 0 2 within Device B in turn 
leads to its higher amperometric response (27-30%) as illustrated in Figure 3. 

Figure 4. (A) Simulated glucose and H 2 0 2 concentration profiles in the multi-layer sensor 
system consisting of 20 nm GO x /PPD as the first layer, and HAs/Fe as the second layer; 
(B) simulated glucose and H 2 0 2 concentration profiles in the multi-layer sensor system 
consisting of 20 nm GO x /PPD as the first layer, and PSS/PDDA as the second layer. 



6.00 



S 5.95 



2 5.90 



u 
u 

§ 5.85 

u 

2 5.80 



(A) 



le - Device A (HAi/Fe 3 )**•- 
„ ,. . ...■>•, 



Glucose Profile - 

H 2 0 2 Profile - Device A (HAs/Fe J 



3.5 
3.0 
2.5 
H 2.0 
1.5 



^ 6.00 



0 20 40 60 80 100 120 

Distance From Electrode Surface (nm) 



3 

e 

o 



1.0 © 
U 

- 0.5 O 
= 



0.0 



5.95 



~ 5.90 



■~> 
s 

O 

u 



5.85 



5.80 



■ Glucose Profile 
Device B (PSS/PDDA) 

• H202 Profile 
Device B (PSS/PDDA).** 



(B) 



» ■ 



3.5 

3.0 5 
2.5 ^ 

_o 

2.0 « 
la 

a 

1.5 H 

u 
c 

-a 

0.5 

s 

o.o 



0 5 10 15 20 25 

Distance from Electrode Surface (nm) 



Figure 5. 90% saturation current response time vs. glucose concentration of Device A and 
B, employing a 20 nm thick GO x /PPD layer. 
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The unwanted build-up of H2O2 within Device B is also expected to increase the sensor response 
time considering the longer times needed to establish equilibrium. This increase in response times adds 
to the already longer response times resulting from the use of tighter (less glucose permeable 
membranes). With this in mind, the response times (tgo%; time taken to reach 90% of the saturation 
current value as obtained from the simulated sensor transient response, Figure S3 of Supplementary 
Information) of the devices have been investigated as a function of glucose concentration and type of 
outer LBL membrane. Figure 5 shows the t9o% response times for the two devices under hypoglycemia 
(2 mM glucose), normal physiological glucose concentrations (6 mM glucose) and hyperglycemia 
(20 mM glucose). For both devices at any glucose concentration, the response time of Device A was 
ca. 35% less than that of Device B, as expected. 

Overall the results of Figures 4 and 5 indicate one has to utilize a tighter PSS/PDDA membrane for 
higher sensitivities (e.g., for glutamate detection which are in the sub-uM levels [35]) even though it 
can lead to longer response times (due to build-up of H2O2) necessitating a means to remove H2O2 
buildup within the sensor. This can be achieved by methods including: (i) increasing the catalytic 
activity of the Pt working electrode [22] (ii) introducing an additional F^C^-consuming layer 
(e.g., Catalase enzyme) [11], and (ii) grading the GO x concentrations in the PPD membrane. Such 
schemes retain the outer boundary conditions of the currently utilized model (Figure 2), however 
the inner boundary conditions must be modified in order to include the associated diffusion and 
reaction-diffusion equations. 

Having established the importance of outwards diffusion of H2O2 on sensor performance, the 
enzyme layer thickness was additionally simulated at 15 nm and 25 nm, in an effort to determine if the 
observed phenomenon is enzyme-thickness dependent. Figure 6 shows the peak H2O2 concentration in 
the GO x /PPD layer of Device A and B at 15 nm (Figure SI of Supplementary Information), 20 nm 
(Figure 4) and 25 nm (Figure S2 of Supplementary Information) thick enzyme layers. Intuitively, it is 
seen that there is increase in peak H2O2 concentration at increasing enzyme thicknesses, due to the 
increase in GO x -induced H2O2 turnover. Moreover, for any enzyme thickness, the H2O2 concentration 
for Device B is ca. 50-70% larger than Device A due to the pronounced buildup of H2O2 in the 
enzyme layer. 

Figure 6. Peak H2O2 concentration of Device A and B in the GO x /PPD layer as a function 
of enzyme thickness (15, 20 and 25 nm). 
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Following the analysis of peak H2O2 concentration as a function of varying GO x /PPD thickness, the 
t9o% response times of these thicknesses was investigated. Figure 7 shows the response times of the 
devices as a function of enzyme thickness. As expected, for any enzyme thickness, the response times 
of Devices B is higher than that of Devices A, due to the higher build-up of H 2 0 2 (shown in Figure 6). 

Figure 7. 90% saturation current response time vs. glucose concentration of Device A and 
B, as a function of GO x /PPD thickness (15, 20 and 25 nm). 
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4. Conclusions 



In this contribution, a multi-layered computational model has been utilized to compare the 
performance of an enzymatic glucose sensor coated with two kinds of LBL assembled outer membranes, 
namely humic acids/ferric cations and poly(styrene sulfonate)/poly(diallyldimethylammonium chloride) 
outer membranes. The diffusion profiles of glucose (S) and H2O2 (P) have been computed under 
specified boundary conditions for sensors coated with high glucose permeable (HAs/Fe 3 ) and low 
glucose permeable (PSS/PDDA) membranes. Sensors coated with PSS/PDDA membranes, although 
have low inwards flux, have displayed high sensitivity due to low outwards diffusion of H2O2 in all 
simulated enzyme thicknesses. This trend in simulation corroborates our previous experimental finding 
of the role H2O2 on sensor sensitivity [10] and re-signifying the need for outer membrane optimization 
not only in terms of glucose diffusion, but also for H 2 0 2 diffusion. Moreover, it is also found that this 
low outer diffusion of H2O2 through the less glucose permeable PSS/PDDA membranes leads to longer 
response times thereby alluding to the necessity for a means to remove H2O2 buildup within an outer 
membrane-coated biosensor. Currently, we are investigating the efficacy of this model in a variety of 
physiological environments as well as understanding the need for gradients, if any. Additionally, even 
though H 2 0 2 serves as the key contributor to the dynamics of these devices, the role of H 2 0 2 in an 
environment of varying oxygen levels is currently under investigation, and will be the topic of future 
publications. 

The study presented here, even though only two LBL membranes were compared, could be 
extended to other designs that utilize a different outer membrane or to designs that employ a series of 
stratified membranes and could be used a guide for optimization of outer layers. 
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